Entanglement and localization of wavefunctions 
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We review recent works that relate entanglement of random vectors to their localization properties. In partic- 
ular, the linear entropy is related by a simple expression to the inverse participation ratio, while next orders of 
the entropy of entanglement contain information about e.g. the multifractal exponents. Numerical simulations 
show that these results can account for the entanglement present in wavefunctions of physical systems. 
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I. INTRODUCTION 

Quantum mechanics has always seemed puzzling since its 
first construction in the first half of the XXth century. Many 
properties are different from the world of classical physics in 
which our intuition is built. The development of quantum in- 
formation science in the last decades has exemplified this as- 
pect. Indeed, it was realized that it is in principle possible to 
exploit the features of quantum mechanics to treat information 
in a different way from what a classical computer would do. 
In this context, the specific properties of quantum mechanics 
are put forward as new resources which enable to treat infor- 
mation in completely new ways. 

One of the most peculiar properties of quantum mechan- 
ics is entanglement, that is the possibility to construct quan- 
tum states of several subsystems that cannot be factorized into 
a product of individual states of each subsystem. Such en- 
tangled states are the most common in quantum mechanics, 
and they display correlations which cannot be seen in a clas- 
sical world, exemplified by e.g. the Einstein-Podolsky-Rosen 
"paradox". Entanglement is also a resource for quantum infor- 
mation (see [1] and references therein), and has been widely 
studied as such in the past few years. 

Despite intensive work, entanglement remains a somewhat 
mysterious property of physical systems. The structure of en- 
tanglement of systems even with small numbers of particles is 
hard to characterize. Even properly measuring the entangle- 
ment present in a system is difficult for mixed states. This is 
all the more important since recent results have shown that (at 
least for pure states) if a process creates a sufficiently low level 
of entanglement, it can be simulated efficiently by a classical 
computer yfl . This gives a limit on the speedup over classical 
computation a quantum computer can achieve, and also gives 
rise to interesting proposals for building classical algorithms 
simulating weakly entangled quantum systems |3fl. 

In this paper, we review recent results we obtained (details 
can be found in 0>H1]), which concern the relationship of en- 
tanglement to localization properties of a quantum state. Our 
strategy is to consider «-qubit systems, and to study entan- 
glement of quantum states relative to their localization prop- 
erties in the 2"-dimensional Hilbert space in the computa- 
tional basis. We obtain analytical results for random states, 
that is ensemble of quantum states sharing some properties. 
Such random states have been recently studied in the liter- 



ature. They are interesting in themselves, since it has been 
shown for example in quantum information that they are use- 
ful in various quantum protocols @|. This motivated a re- 
cent activity in the quantum information community to try 
and produce efficiently such random vectors or random op- 
erators through quantum algorithms 07|], and to characterize 
their entanglement properties [8]. In addition to their intrin- 
sic usefulness, random states are important since they can de- 
scribe typical states of a "complex" system. For example, it 
has been known for some time now that random vectors built 
from Random Matrix Theory (RMT) can describe faithfully 
the properties of quantum Hamiltonian systems whose clas- 
sical limit is chaotic, and more generally of many complex 
quantum systems Such random vectors are ergodic, and 
the entanglement they contain has been calculated some time 
ago II 1 OL 1 1 lfl . However, in many quantum systems, the wave- 
functions are not ergodic but localized. This can correspond 
to electrons in a disordered potential, which are exponentially 
localized due to Anderson localization. It can also be seen 
in many-body interacting systems, where the presence of a 
moderate interaction can lead to states partially localized in 
energy. Some systems are in a well-defined sense neither er- 
godic neither localized: they correspond to e.g. states at the 
Anderson transition between localized and delocalized states, 
and can show multifractal properties lfl2ll . 

In this paper, we calculate the amount of entanglement 
present in ensembles of random vectors displaying these var- 
ious degrees of localization. Besides generalizing the result 
for RMT-type random vectors, this gives the entanglement 
present in a "typical state" of such localized or partially local- 
ized systems. This enables to estimate the complexity of sim- 
ulating such systems on classical computers, but also sheds 
light on the entanglement itself, since in these cases it is re- 
lated through simple formulas to quantities characterizing the 
degree of localization of the system. 

Our results show that for random vectors which are local- 
ized on the computational basis, the linear entropy which ap- 
proximates the amount of entanglement in the vector is sim- 
ply related to the Inverse Participation Ratio (IPR), a popular 
measure of localization. The next term in the approximation is 
related to higher moments, and in particular to the multifractal 
exponents for multifractal systems. In order to assess the use- 
fulness of these results to physical systems, we compare them 
to the entanglement numerically computed for several models. 
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After a general discussion on entanglement of random vectors 
(section II), we consider the entanglement of one qubit with 
the others (section III), and give explicitly the first and second 
order of the expansion of the entropy of entanglement around 
its maximum. Section IV generalizes these results to other bi- 
partitions, and section V compares the formula obtained with 
the numerical results for two physical systems. Section VI 
considers the physically important case of vectors localized 
not on a random subset of the basis vectors, but on a subset 
composed of adjacent basis vectors (that is the states are lo- 
calized on computational basis states which are adjacent when 
the basis vectors are ordered according to the number which 
labels them), showing that the results become profoundly dif- 
ferent. Section VII presents the conclusions. 



II. ENTANGLEMENT OF RANDOM VECTORS 

Random vectors are ensembles of vectors whose compo- 
nents are distributed according to some probability distribu- 
tion. If for example the system considered is composed of 
n qubits, the Hilbert space is of dimension N = 2". If the 
two states of a qubit are denoted |0) and |1), each state in 
the computational basis corresponds to a sequence of and 1 
and thus can be labelled naturally by a number between and 
2" - 1, and quantum states can be expanded as \ifr) = 
Random vectors distributed according to the uniform measure 
on the Af-dimensional sphere describe typical quantum states 
of the n qubits. Such states are ergodically distributed in the 
computational basis, and their entanglement has already been 
studied in ifTol [Till . In this paper, we are interested in ran- 
dom vectors which are not ergodically distributed. Ensembles 
of such states will be characterized by localization properties. 
The simplest example of such localized random vectors can be 
constructed by taking M components (M < N) with equal am- 
plitudes and uniformly distributed random phases, and setting 
all the others to zero. The random vectors will all be exactly 
localized on M basis states. A more physically relevant ex- 
ample consists in still choosing M < N nonzero components, 
and giving them the distribution of column vectors ofMxM 
random unitary matrices drawn from the Circular Unitary En- 
semble of random matrices (CUE vectors). In general our re- 
sult will be averaged both over the distribution of the nonzero 
components and the position of these nonzero components in 
the computational basis. This corresponds to classes of ran- 
dom vectors sharing the same localization length. Our results 
will in fact generalize to any such distribution of random vec- 
tors whose localization properties are fixed. In addition, we 
shall see that if we impose that the distribution of the indices 
i of nonzero components iff, is such that they are always adja- 
cent in the computational basis (i.e. the indices i are consecu- 
tive integers), the results change drastically. 

The localization properties of the random vectors can be 
probed using the moments of the distribution 



Pq 



Z i* 1 



2<l 



(D 



The second moment is j?2 = 1 /£, where £ is the Inverse Par- 
ticipation Ratio (IPR) which is often used in the mesoscopic 
physics literature to measure the localization length. Indeed, 
for a state uniformly spread on exactly M basis vectors, one 
has £ = M. The scaling of p2 and higher moments with the 
size also probes the multifractal properties of the wavefunc- 
tion. 

The random states we consider are built on the N- 
dimensional Hilbert space of a n-qubit system with N = 2". 
We are interested in bipartite entanglement between subsys- 
tems defined by different partitions of the n qubits into two 
sets. In general, bipartite entanglement of a pure state belong- 
ing to a Hilbert space "Ha ® He is measured through the en- 
tropy of entanglement, which has been shown to be a unique 
entanglement measure 111 311 . We consider pure states belong- 
ing to 'Ha ® 'Hb where "Ha is a set of v qubits and fig a set of 
n — v qubits. If pa = trg|^)(^| is the density matrix obtained 
by tracing out subsystem B, then the entropy of entanglement 
of the state if/ with respect to the bipartition (A, B) is the von 
Neumann entropy of pA, that is S — -tr(p^ log 2 pa)- 



III. ENTANGLEMENT OF ONE QUBIT WITH ALL THE 
OTHERS 

To obtain an approximation for the entropy, one can expand 
S around its maximal value. In the case of the partition of the 
n qubits into 1 and n — 1 qubits, the entropy can be written as 
a function of r, with 



t = 4detpA 



(2) 



(in the case of 2 qubits this quantity is called the tangle 
and corresponds to the square of the generalized concurrence 
01). One has 



5(t) = h 



(\+ VT^ 



(3) 



where h(x) = — x log 2 x - (1 — x) log 2 ( 1 - x). The series expan- 
sion of S (t) up to order m in (1 - r) reads 
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The first order corresponds to r itself up to constants and 
its average over the choice of the (1, n - 1) partition is known 
as the linear entropy or Meyer- Wallach entanglement Q lfl5ll . 
Our results show that for our class of random vectors, the av- 
erage linear entropy is given by 
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This formula was obtained first by considering a random 
vector which is nonzero only on M basis vectors among N, 
and summing explicitly the combinatorial terms. It can also 
be obtained in a more general setting by taking M — N and 
summing up all the localization properties of the vector in the 
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IPR £ alone. For any (l,n - 1) partition of the n qubits, the 
components of the vector can be divided in two sets according 
to the value of the first qubit. Assuming no correlation among 
these sets enables to get Eq. (0 (details on the calculations 
can be found in IJ]). 

It is interesting to compare this formula with a similar 
one obtained in iflol using different assumptions, in partic- 
ular without average over random phases. The formula ob- 
tained relates entanglement to the mean inverse participation 
ratio calculated in three different bases, a quantity that is more 
general but often delicate to evaluate. In our case, the addi- 
tional assumption of random phases enables to obtain a for- 
mula which involves only the IPR in one basis, a quantity that 
can be easily evaluated in many cases. For example, it en- 
ables to compute readily the entanglement for localized CUE 
vectors. However there are instances of systems (e.g. spin sys- 
tems) where these different formulas give the same results. 

In particular, our formula (0 allows to compute (t) e. g. 
for a CUE vector localized on M basis vectors; in this case 
£ = (M + l)/2, and we get 



<r> = 



M-lN-2 
M+lN-l' 



(6) 



In [10], (t) was calculated for non-localized CUE vectors of 
length N, giving (r) = (N - 2)/(N + 1). Consistently, our 
formula yields the same result if we take M — N. For a vec- 
tor with constant amplitudes and random phases on M basis 
vectors, £ = M and 



<T> = 



M-lN-2 
M N-l' 



(7) 



Thus the first order of the expansion, which gives the main 
features of the entanglement, has very simple expressions in 
term of system parameters. 

The next order in the expansion can be obtained by sim- 
ilar methods that we do not detail here (see [5] for details); 
summing up all terms involved in r 2 we get 

<t 2 > = N(N-2)(N 2 -6N+ 16)cim (8) 
+4N(N - 2)(N - 4)c 2 i i + 4V(V - 2)c 22 



with 
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This gives the next order of the entropy of entanglement 
in terms of the moments up to order 4 of the vector. What 
this means is that at this order, the average entanglement of 
random vectors with fixed moments will be related to them 
through ||8). Although more complicated than ([5j, the for- 
mula indicates that e.g. for states having multifractal proper- 
ties, since moments scale with system size according to quan- 
tities called multifractal exponents, the behavior of the entan- 
glement at this order will be also controlled by these multi- 
fractal exponents. 



The nth order of the expansion dU can similarly be obtained 
and has been derived in 051]. It is interesting to note that in the 
case of a CUE random vector of size N, resummation of the 
whole series for S (r) yields, after some algebra, 



(S(t)) = 
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which has been obtained earlier by a different method [11]. 

A general conclusion obtained from these formulas is that 
the entanglement associated to such bipartition (1, n - 1) goes 
to the maximal value for large N and large even if £ grows 
more slowly than N. For fixed it tends for large N to a 
constant nonzero value which depends on We will see in 
Section VI that this result can change drastically if we impose 
a localization on fixed locations in Hilbert space. 



IV. ENTANGLEMENT OF RANDOM VECTORS: OTHER 
PARTITIONS 

Up to now we have considered the entanglement of one 
qubit with all the others, i.e. the (I, n-l) partition of n qubits. 
What about bipartite entanglement relative to other biparti- 
tions (v, n - v), where v is any number between 1 and n — 1? 
In this case, it is convenient to define the linear entropy as 
S l = grjCl _ a "Pi)' wnere d = dim^A < dim 'Kb. The scal- 
ing factor is such that S l varies in [0,1]. 

A similar calculation as above enables then to obtain the 
first order of the mean von Neumann entropy. It is given by 
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with p2 = 1/f, which generalizes Eq. ©. 

Higher-order terms can be obtained as well, although the 
calculations become tedious. To this end, the entropy S = 
-tr{pA log 2 Pa) is expanded around the maximally mixed state 
po = l/2\ as 

S = V+ ri- f ^ZTMiPA -P(»)" +1 ), (12) 

In 2 n(n + 1) 



and the traces can be evaluated as sums over correlators of 
higher moments |5J] . 

We remark that again the linear entropy (TTTb tends to the 
maximal possible value when N and £ become large, as for 
the (l,7i - 1) partition. 



V. ENTANGLEMENT OF RANDOM VECTORS : 
APPLICATION TO PHYSICAL SYSTEMS 

In order to test these results on physical systems, we com- 
pared them to numerical results obtained from different mod- 
els. 

The first one corresponds to a diagonal Hamiltonian matrix 
to which a two-body interaction is added. 
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FIG. 1: Scaled mean linear entropy (t)(N - 2)/N of U3\ vs mean 
IPR for 5 = Ao, n = 10 (blue circles) and n = 11 (green squares). 
Red line is the theory, stars the data for n = 10 with random shuffling 
of components (from jj]). 
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FIG. 2: Mean entropy of entanglement 5 for different bipartitions 
(v, n — v) as a function of the mean IPR. Left: eigenvectors of d 14b 
with y = 1/3; the average is taken over 10 6 eigenvectors. Right: 
eigenvectors of J 1 3 fc with 6 = Ao and J/6 = 1.5; average over « 
3 x 10 5 vectors. Triangles correspond to v = 1, squares to v = 2 and 
circles to v = n/2, with n = 4— 10. Black symbols are the theoretical 
predictions for the mean value of S (obtained from Eq. I ll It and green 
(gray) symbols are the computed mean values of the von Neumann 
entropy (from [5J|). 
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This system can describe a quantum computer in presence 
of static disorder ifnll . Here the <x,- are the Pauli matrices for 
the qubit i, the energy spacing between the two states of qubit 
i is given by T,, which are randomly and uniformly distributed 
in the interval [Ao - 6/2, Ao + 6/2], and Jjj uniformly dis- 
tributed in the interval [-J, J] represent a random static inter- 
action. Entanglement of eigenvectors of this Hamiltonian was 
already considered in a different context in 111 811 . It is known 
lfl7ll that in this model a transition to quantum chaos takes 
place for sufficiently large coupling strength J. In this regime, 
eigenvectors of dTJb are spread over all noninteracting eigen- 
states (those of ( fT3l ) for J — 0, which coincide to the com- 
putational basis), but in a certain window of energy, and are 
distributed according to the Breit-Wigner (Lorentzian) distri- 
bution. Thus these wavefunctions are distributed among a cer- 
tain subset of the computational basis, although they are not 
strictly zero outside it, and the distribution is not uniform, but 
rather Lorentzian. Nevertheless, our data show (see Figs. 1,2) 
that the behavior of the bipartite entanglement of eigenvec- 
tors of this model is well described by the results (O and ( fTTT i 
derived for random vectors. The agreement becomes very ac- 
curate if the eigenvector components are randomly shuffled to 
lower correlations. 

We also considered another model, based on NxN matrices 
of the form 



1 



a 2mNy 



N 1 



( ,2iMk-l+Ny)/N ' 



(14) 



where <f>k are random variables independent and uniformly dis- 
tributed in [0, 2n[ and y is a fixed parameter. This model intro- 
duced in fl9ll is the randomized version of a simple quantum 
map introduced in 112011 . The eigenvectors of (fl4l have multi- 
fractal properties Il2lll for rational y. The results of Figs. 2,3 
show that again the results for random vectors describes very 
well the entanglement for this system for randomly shuffled 
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FIG. 3: Relative difference of the entropy of entanglement ([3} and its 
successive approximations S „, (m = 1, 2) with respect to the number 
of qubits for eigenvectors of j 14b for (left) 7 = 1/3 and (right) 7 = 
1/7. The average is taken over 10 7 eigenvectors, yielding an accuracy 
< 10" 6 on the computed mean values. Green triangles correspond 
to the first order expansion 5 1 , blue squares and red circles to the 
second order expansion 52- The difference between the latter two is 
that for blue squares (pi) appearing in Eq. $9$ has been replaced by 
{pi) 2 yielding a less accurate approximation. Dashed line is a linear 



fit yielding 1 - <5 1)/<5> 
(from d). 



Af- 



ter 7 = 1/3 and A^ 1 - 58 for 7 = 1 /7 



components, and that even the first order is already a good 
approximation. 



VI. ENTANGLEMENT OF ADJACENT RANDOM 
VECTORS 



In the preceding sections we discussed formulas for entan- 
glement of ensembles of random vectors where the compo- 
nents over each basis vector are independent. If we relax this 
assumption, the result may change. An important particular 
case corresponds e.g. to random vectors localized on M com- 
putational basis states which are adjacent when the basis vec- 
tors are ordered according to the number which labels them 
(again, if the two states of a qubit are denoted |0) and |1), each 
state in the computational basis corresponds to a sequence of 
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FIG. 4: Mean linear entropy (t) of partitions (l,n - 1) vs number of 
qubits for the one-dimensional Anderson model with disorder from 
top to bottom w = 0.2 (blue), 0.5 (red), 1.0 (green), 1.5 (magenta), 
2.0 (cyan), and 2.5 (orange). Average is over 10000 eigenstates. 
Solid lines are the C/n fits of the tails (from @]). 



is a diagonal operator whose elements e, are Gaussian ran- 
dom variables with variance w 2 , and V is a tridiagonal matrix 
with non-zero elements only on the first diagonals, equal to 
the coupling strength, set to 1 . It is known that eigenstates of 
this system, which modelizes electrons in a disordered poten- 
tial, have envelopes of the form exp(-|x - Xo\/l), where / is 
the localization length. It was shown in [ 22] that this model 
can be simulated efficiently on a quantum computer, and the 
wavefunction of the computer during the algorithm will be 
localized on adjacent basis vectors, which correspond to the 
position of vertices. Figure 4 shows that the asymptotic be- 
havior of the linear entropy of the eigenstates (with all corre- 
lations left between components, i.e. no random shuffling) is 
well captured by Eq. ( fl~5b . 

Thus random vectors localized on adjacent basis vectors 
correspond to a drastically different behavior compared to the 
vectors of section III : indeed, for fixed £ the entanglement (at 
least the linear entropy) always tends to zero for large N, even 
if it does it rather slowly (as ~ 1/ lnAO- 



and 1 and thus can be labelled naturally by a number be- 
tween and 2" - 1). In this case, we had to use combinatorial 
methods; summing all contributions together we get for the 
linear entropy of ( 1 , n - 1 ) partitions 
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M - 2 2{2 r « 
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where ro is such that 2 r ° 

v-2 2 „(Jl 



< M <2 r ° and^ r (*) = ;tv(2 r+1 - 
x) = x l - \x{x 2 - l)/2 r for < x < 2 r . Equation ([151) is an 
exact formula for M < N/2. For fixed M and n — > oo, n(Q) 
converges to a constant C which is a function of M and For 
M = 2'°, ro < n, Eq. ( TBI ) simplifies to 



<T> = 



f(ro + f)M 2 -2(r -l)M-f 
M( M - 1) 

^)K>)];- 



(16) 



Numerically, this expression with ro = log 2 M gives a very 
good approximation to Eq. (TT~5T > for all M. 

Equation ([T5T l is exact for e.g. uniform and CUE vectors, 
and can be applied even if the vector is not strictly zero out- 
side a M-dimensional subspace. Indeed, for A^-dimensional 
CUE vectors with exponential envelope exp(-x/Z), (Q) is in 
excellent agreement with Eq. (TT3T > with £ = I and M = 2%. 

In order to compare these findings to those of a physical 
system with such a property of localization on adjacent ba- 
sis vectors, Fig. 4 shows the theory Eq. (TT3T > together with 
the entropy for the one-dimensional Anderson model. This 
model corresponds to a one-dimensional chain of vertices 
with nearest-neighbor coupling and randomly distributed on- 
site disorder, described by the Hamiltonian Hq + V. Here Hq 



VII. CONCLUSION 

The results above indicate that the entanglement properties 
of random vectors can be directly related to the fact that they 
are localized, multifractal or extended. The numerical simu- 
lations for different physical systems show that these results 
obtained for random vectors describe qualitatively the entan- 
glement present in several physical systems, and reproduce it 
accurately if correlations between components of the vector 
are averaged out. 

Thus the results are interesting to predict the amount of en- 
tanglement present in random vectors, and also can be applied 
to physical systems for which such random vectors describe 
typical states. This gives insight on the difficulty to simu- 
late classically such systems, since systems with low amounts 
of entanglement can be simulated classically efficiently. This 
also can be applied to estimate the changes in entanglement at 
a quantum phase transition l23ll . in particular for the Ander- 
son transition between localized and extended states (see JH 
for more details). Additionally, this gives also insight on the 
nature of entanglement itself by relating it to simple physical 
properties of the system. 
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